QUANTUM HAMILTON- JACOBI THEORY 



=3 . 



Marco Roncadelli 

INFN, Sezione di Pavia, Via A. Bassi 6, 1-27100 Pavia, Italy, 
and Dipartimento di Fisica Nucleare e Teorica, Universitd di Pavia, Ital^ 

L. S. Schulman 



l^- , rnysics uepartment, uiartison university, 

Potsdam, New York 13699-5820, UsM 
(N 1 (Dated: February 2, 2008) 



Quantum canonical transformations have attracted interest since the beginning 



Q | of quantum theory. Based on their classical analogues, one would expect them 

to provide a powerful quantum tool. However, the difficulty of solving a nonlinear 
operator partial differential equation such as the quantum Hamilton-Jacobi equation 
(QHJE) has hindered progress along this otherwise promising avenue. We overcome 
Qh| this difficulty We show that solutions to the QHJE can be constructed by a simple 

prescription starting from the propagator of the associated Schrodinger equation. 
Our result opens the possibility of practical use of quantum Hamilton-Jacobi theory. 



As an application we develop a surprising relation between operator ordering and 
the density of paths around a semiclassical trajectory. 

PACS numbers: 03.65.Ca 
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Canonical transformations play a central role in classical mechanics jlj . From the earliest 
days of quantum mechanics, the importance of quantum canonical transformations (QCT's) 
has been recognized and their properties have been systematically investigated by Jor- 
dan [3], London Dirac 0] and Schwinger among others. Schwinger's framework — 
based on the Quantum Action Principle [7| — provides the most suitable context to define a 
QCT as q — > Q,p — > P, H(q, p, t) — > K(Q, P, t), where all canonical variables pertain to the 
same dynamical system S with degrees of freedom 



X 

Owing to the formal similarities between classical and quantum mechanics, QCT's closely 
resemble their classical counterparts. In particular, one of four possible sets of independent 
canonical variables (q,Q), (q,P), (Q,P), (P,P) must be selected to represent a QCT ex- 
plicitly, and we denote by W(q,Q,t), W(q,P,t), etc., the associated operator generating 
functions. Choosing the set (g, Q), a QCT can be written as (1 < i < N) 

Pi = £rW(q,Q,t) , (1) 
dqi 

() W(q,Q,t) , (2) 



dQi 
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K{Q, P, t) = H(q, p, t) + ^W{q, Q, t) . (3) 

The presence of noncommuting operators in the generating function makes QCT's differ- 
ent from the classical ones and is ultimately responsible for the difference between classical 
and quantum mechanics (||. As emphasized by Jordan and Dirac, the resulting operator- 
order ambiguity should be fixed by enforcing well-ordering: operators represented by capital 
letters should always stay to the right of those labelled by lower case letters. This means 
that W(q,Q,t) should have the structure 

W(q,Q,t)=Y,fa(q,t)g a (Q,t) , (4) 

a 

for suitable functions /«(•) and g a (-). Throughout, we will suppose that operator generating 
functions are well-ordered. Note that with well- ordering a quantum generating function like 
W(q,Q, t) is uniquely defined by the replacements q — > q, Q — ► Q in a given c-number 
function W(q,Q,t) [fl. 

As in classical mechanics, the quantum time evolution is described by a canonical trans- 
formation bringing the canonical variables in the Heisenberg picture q(t), p(t) to constant 
values at some initial time to- In addition, q(t), p(t) can be derived from Eqs. ([I]) and (j2J), 
provided that the transformed Hamiltonian vanishes. As a consequence, the operator gen- 
erating function W(q,Q,t) obeys the operator quantum Hamilton- Jacobi equation (QHJE) 

H U ^W(q, Q, t),tj+ ^W(q, Q,t) = 0. (5) 

W(q, Q,t) should be a complete solution of Eq. (jSJ), i.e., it should depend on iV indepen- 
dent "integration constants" Qi. As in classical mechanics, the operator Hamilton- Jacobi 
equation, Eq. (J5J), provides an independent formulation of the theory. Yet the formidable 
difficulty of finding solutions to this nonlinear operator partial differential equation has 
hindered progress along this otherwise promising avenue. 

Our aim is to show that this stumbling block can be sidestepped, thereby opening the 
way to exploiting the operator QHJE as a calculational tool. As we will demonstrate, 
the solutions to the operator QHJE arise by a simple prescription from the solutions of 
the Schrodinger equation for the same Hamiltonian. In particular, the operator generating 
function W(q,Q,t) arises from the quantum propagator. Implications of our result will be 
discussed after we have completed its demonstration. 

We are concerned throughout with the general Weyl- ordered Hamiltonian [ll| 

H(q,p,t) = -dijifipiPj + Piaijifipj + ^PiPjdijiq) 
+bi(q)pi + Pih{q) + c(q) , 

(6) 

where dij(-), bi(-), and c(-) are functions of and summation over repeated Latin indices 
for the degrees of freedom of S is understood. Employing the shorthand W = W(q,Q,t), 
Eq. USD reads 

1 ,^dWdW dW ,^dW ldWdW , A , , ,^dW 
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Since we are looking for the relationship between the operator QHJE and the Schrodinger 
equation, we turn Eq. (JTj) into a c-number partial differential equation. Hence we sandwich 
Eq. ((Tj) between (q\ and \Q), finding 



~-aij(q)(q 



1, 
+ 2 {q 



dW dW 



dqi dqj 
dW dW 



Q) + (q 



dW , ^ dW 



Q) + bi(q)(q 



dqj 
dW 



+(q 



dW 



bi{q) 



Q) + c(q)(q\Q) + (q 



dqi 
dW 



dt 



Q) 
Q) 

Q) =o 



(8) 



To evaluate the matrix elements in Eq. (jSJ), we make repeated use of the canonical commu- 
tation relations. In this connection, we recall that for an arbitrary function G(-) 



[G(q),Pi] = ih 



9G(q) 
dqi 



(9) 



By inserting Eq. ((T]) into Qj, we obtain 



dW r(-\ r(^ dW -K dG (fi 

-w, G{q)=G{q) wr lh ^q- 



(10) 



We begin by taking G(q) = bi(q), G(q) = Oy(g) and G(q) = da,ij(q)/dqj. Accordingly, Eq. 
(|T0|) allows us to rewrite Eq. (jSj) as 



2a ij (q)(q 



dW dW 



dqi dqj 



Q) + 2[b i (q)-ih^^-)(q 



dW 



+ c (q) - 



..dbi(q) fr 2 d 2 a ij (q) 



dqi 



2 dqidq. 



(q\Q) + (q 



dqi 
dW 



dt 



Q) 

Q) = 0. 



(11) 



At this point, we denote by W(q,Q,t) the c-number function that uniquely produces 
W(q, Q, t) by the substitution q — > q, Q — > Q. Explicit use of Eq. (0J yields 



(q\W\Q) = W(q,Q,t)(q\Q), 
dW 



dt 
dW 



and furthermore 



dW dW 



dqi dqj 



dW(q,Q,t) , 

Q) = — -Q t — \q\Q) 

8W{q,Q,t) , 
Q) = dgi (q\Q) 



(12) 
(13) 

(14) 



dqi 



t) 



Q) 



(15) 
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What remains to be done is to disentangle Eq. (|T5|) . To this end, we first take G(q) = 
dfp{q,t)/dqj in Eq. ([10]) to get 

\- df a (q,t) (A dfpjq^t) _ df p {q,t) ^ df a {q,t) ~ 
^ dq, 9a[Li ' ] dq, ~ dq 3 ^ dq, 9aW) 

^d 2 f (q,t) , , 

' (16) 

We next multiply Eq. f|T6|) by gp(Q,t) on the right and sum over a, thereby obtaining 

E %M 9 „ ( Q, t) ^) 9 ,(Q, t) 



o,/3 



dfp{q, t) df a (q, t) A A ^ d 2 fp(q,t) A 

f^—9a(Q,t) g/3 (Q,t) -ih^ g p {Q,t) , (17) 



a,/3 1J 11 /3 

which allows us to rewrite Eq. (1T51) as 
, OWdW 



dq-i dq~j 

As a consequence, Eq. ffTTl) takes the form 



2a . -(g) ( dW M,t) dW(q,Q,t) _ ^Wjq^t) 
tj V % dqj dqidqj 

+2 ( bi (q) - in^M") + c{q) _ JbM) 



dqj J dqi ' % 

V cPaM dW(q,Q,t) _ n 
2 %c% + <9t 1 } 

This derivation makes it natural to regard Eq. (1T91) as the c-number QHJE associated with 
the operator QHJE (JSJ) for 5 described by the quantum Hamiltonian 
The physical significance of Eq. (|T9|) becomes clear by setting 

Q, t) = exp{(i/h)W(q, Q,t)}. (20) 

A straightforward (if tedious) calculation shows that i[)(q,Q,t) obeys precisely the 
Schrodinger equation associated with the quantum Hamiltonian ([SD in the variables q, t [l2| • 
Hence - thanks to Eqs. (1T21) and (1201) - starting from a solution W(q, Q, t) we get a solution 
ijj(q, Q, t) of the corresponding Schrodinger equation depending on N independent constants 
Qi. We stress that this result holds true even for solutions W(q, t) of the operator QHJE 
that are independent of Q, since all equations from (j7j) onward could have been multiplied 
by f dQ(f>(Q), with 4>(Q) arbitrary. What is more important for us, the argument can be 
turned around, because W(q,Q,t) can be uniquely obtained from W(q,Q,t) by enforcing 
well- ordering. Therefore, from a solution ip(q,Q,t) of the Schrodinger equation depending 
on N independent constants Qi we get W(q, Q, t), and from any particular solution ip(q,t) 
we can construct a particular solution W(q, t). 
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So far, we have focused on showing that W(q, Q, t) satisfies a certain differential equation. 
As we demonstrate below, by use of appropriate boundary conditions we get more specific 
information. Namely, the solution of Schrodinger's equation that results from the operator 
generating function W(q,Q,t) is precisely the quantum propagator K(q,Q,t) [13j . Since 
any solution of the Schrodinger equation arises by convolving an arbitrary wave function 
with the propagator, we conclude that any solution of the operator QHJE can ultimately 
be constructed in terms of the propagator. 

This will allow solutions of and approximations to the operator QHJE to be obtained, 
since a wealth of information is available on the corresponding solutions to Schrodinger's 
equation. In particular, once an exact or approximate W has been constructed, one can 
obtain the time dependence of operators, using Eqs. §I§ and (j2]). 

We proceed to prove that ip(q,Q,t) = K(q,Q,t). Since both quantities satisfy the same 
Schrodinger equation, which is first order in time, all we need show is that they have the 
same boundary conditions at t — 0. The propagator of course is 5(q — Q) at t — 0. To 
show that ip(q, Q, t) shares this property, we must look at the behavior of W for t — > 0. We 
expect W to generate the identity transformation in the limit t — > 0, but there is a slight 
complication: As in classical mechanics the identity transformation using the (q,Q) 
variables does not have a simple form. 

A way out of this difficulty relies on the observation that for a nonsingular potential 
the solution to the classical Hamilton- Jacobi equation approaches that of the free particle 
for t — > 0. Thus, for sufficiently small t the classical generating function has the form 
F{<1) Q) — m (Q — <?) 2 /2t- We use this to guess the limit of the operator W for t — > 0, and 
from that to obtain the corresponding limit of the c-number function W(q,Q,t). The first 
observation is that as a candidate for the small-t limit of W, the well-ordered operator form 
of F(q, Q) (which contains —2qQ) does not work, which is to say, it does not satisfy Eq. 
To see this in detail — and to see the cure — we assume the following small-t limiting form for 
W 

W=^ t [Q 2 -2qQ + q 2 )+g(t). (21) 

Substituting into Eq. the squaring of W generates a term — (qQ + Qq), rather than 
—2qQ, so that satisfying Eq. §5§ requires 

For small t, one can again neglect the influence of the potential terms and the commutator 
can immediately be deduced from the relation q = Q + Pt/m, the solution of the free particle 
Heisenberg equations of motion. Eq. (1221) now becomes dg/dt = ifijlt and we obtain 

W = -(Q 2 -2qQ + q 2 )+-\nt, for *- , (23) 
-exp(--(Q 2 -2gQ-g 2 )J , for t -> , (24) 

with the constant in Eq. (T2"4"|) arising from integrating dg/dt = ih/2t. It is remarkable that, 
aside from the constant (which is not fixed by W), the gQ-commutation relation has given 
us precisely the correct time-dependence of the propagator. This completes our proof. 

The c-number QHJE f lT9|) has repeatedly attracted interest. For instance, Eq. (|T9l) has 
been derived from a diffeomorphic covariance principle based partly on an SL(2, C) algebraic 
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symmetry of a Legendre transform [14j]. Alternatively, Eq. (fl9l) has been taken as the 
starting point of a classical-like strategy to define c-number quantum action-angle variables 
in quantum mechanics [15|. We also remark that a variant of Eq. (fl9l) has been derived 
within the phase-space path- integral approach to quantum mechanics [161 ]. 

We next show the power of the relation we have just developed between W(q,Q,t) and 
K(q,Q,t), using, as suggested above, known information about the propagator. Consider 
a situation where the semiclassical approximation is valid and there is but one classical 
path between the initial and final points. Then in this approximation, as is well-known, 
K(q,Q,t) = const • ^/det d 2 S/dqdQ exp(iS(q, Q,t)/h), with S(q,Q,t) Hamilton's principal 
function (a solution of the classical Hamilton- Jacobi equation). It then follows from our 
result that W(q,Q,t)\ wo = S(q, Q,t)\ wo — |iMogdet d 2 S/dqdQ\ wo , where 'WO" stands for 
"well-ordered." Now imagine that this expression is inserted in Eq. (J5]). If not for the well 
ordering, S alone would solve the equation. Therefore we conclude that the effect of the well- 
ordering is precisely to demand the presence of the additional term, |iMogdet d 2 S/dqdQ 
(where 'WO" has been dropped because there is already an h in the expression). But that 
additional term (famously) has a meaning of its own: it goes back to van Vleck and represents 
the density of paths along the classical path; it plays an essential role, for example, in the 
Gutzwiller trace formula. What our result says is that this density of paths can be thought 
of as arising from the commutation operations necessary to bring S to well-ordered form. 
Thus the purely quantum issue of commuting operators produces a quantity that one would 
have thought is exclusively derivable from classical mechanics. 

We remark that this relation took us completely by surprise. To check it, we worked the 
simplest non-trivial example we could (our proof above comparing the boundary conditions 
for K and W already showed it to be true for the free particle case). Let H = p 2 /2 + V with 
V = Vo©(a/2 — \x\) and x in one dimension. To lowest order in V the action is S(x,y,t) = 
(x—y) 2 /2t—V at/(x—y) for y < —a/2 and x > a/2. We checked our relation, with x — > q and 

y — > Q and with well-ordering implemented by l/(q — Q) = J n °° du exp(— uq) exp(uQ). 

L J WO 

Using the Baker- Campbell-Hausdorff formula and other techniques and keeping only lowest 
order in V and h, indeed the relation checked out! 

In conclusion, we have shown how to construct solutions to the operator quantum QHJE 
starting from the quantum propagator K(q, Q, t) for the same Hamiltonian. Explicitly, once 
K(q, Q, t) is known we get its "complex phase" W(q, Q, t) via Eq. (1201) . Then, by demanding 
well- ordering, the replacement q — > q, Q — > Q uniquely produces the operator W(q,Q,t). 
Alternatively, by convolving K(q, Q, t) with an arbitrary 4>{Q) we produce any solution of the 
Schrodinger equation. Finally, by replacing q — > q in its "complex phase" we get a solution 
W(q,t) of the operator QHJE. While this is obviously true for exact propagators, it also 
enables one to find approximate solutions to the operator QHJE by exploiting approximate 
propagators. In particular we used the semiclassical approximation to the propagator to 
show that the commutation operations establishing well-ordering provide just what is needed 
to get the density of paths around the classical path. This density of paths satisfies a 
continuity equation which, as O'Raifeartaigh and Wipf [13] emphasize, is in a sense of order 
H (even though it involves classical quantities only and has no K in it!). Although our proof 
establishes this surprising relation, there remains the provocative question of understanding 
its intuitive basis. 
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